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The unfolded protein response (UPR) responds to disruption of endoplasmic reticulum 
(ER) function by initiating signaling cascades that ultimately culminate in extensive 
transcriptional regulation. Classically, this regulation includes genes encoding ER 
chaperones, ER-associated degradation factors, and others involved in secretory protein 
folding and processing, and is carried out by the transcriptional activators that are produced 
as a consequence of UPR activation. However, up to half of the mRNAs regulated by ER 
stress are downregulated rather than upregulated, and the mechanisms linking ER stress 
and UPR activation to mRNA suppression are poorly understood. To begin to address 
this issue, we used a "bottom-up" approach to study the metabolic gene regulatory 
network controlled by the UPR in the liver, because ER stress in the liver leads to lipid 
accumulation, and fatty liver disease is the most common liver disease in the western 
world. qRT-PCR profiling of mouse liver mRNAs during ER stress revealed that suppression 
of the transcriptional regulators C/EBPa, PPARct, and PGC-1a preceded lipid accumulation, 
and was then followed by suppression of mRNAs encoding key enzymes involved in fatty 
acid oxidation and lipoprotein biogenesis and transport. Mice lacking the ER stress sensor 
ATF6ct, which experience persistent ER stress and profound lipid accumulation during 
challenge, were then used as the basis for a functional genomics approach that allowed 
genes to be grouped into distinct expression profiles. This clustering predicted that ER 
stress would suppress the activity of the metabolic transcriptional regulator HNF4a — a 
finding subsequently confirmed by chromatin immunopreciptation at the Cebpa and Pgcla 
promoters. Our results establish a framework for hepatic gene regulation during ER stress 
and suggest that HNF4a occupies the apex of that framework. They also provide a unique 
resource for the community to further explore the temporal regulation of gene expression 
during ER stress in vivo. 
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INTRODUCTION 

Originally identified as a program for improving protein fold- 
ing during ER stress, the vertebrate UPR is composed of 
three separate but functionally overlapping pathways that cul- 
minate in transcriptional upregulation of genes involved in 
ER protein folding and processing (Ron and Walter, 2007). 
The ER-resident endoribonuclease IRE1, which is conserved 
among eukaryotes and exists as a [GenBank:NM_023913.2] 
and p [GenBank:NM_012016.2] paralogs in mammals, cat- 
alyzes splicing of Xbpl mRNA [GenBank:NM_013842.3 and 
NM_001271730.1] to remove a 26 base intron and allow 
for the translation of a transcriptional activator of the 
bZIP family. The PERK kinase [GenBank:NM_010121.2] is 



Abbreviations: bZIP, basic leucine zipper; ER, endoplasmic reticulum; TM, tuni- 
camycin; UPR, unfolded protein response; VLDL, very low density lipoprotein. 



metazoan-specific, and it phosphorylates the translation initia- 
tion factor eIF2ct [GenBank:NM_001005509.2] when activated, 
resulting in transient inhibition of protein synthesis but also 
specific translation of Atf4 mRNA [GenBank:NM_009716.2] to 
produce the bZIP transcriptional activator ATF4. Other ER stress- 
independent eIF2a kinases exist, and phophorylation of eIF2oi 
and the attendant consequences of that event are known as the 
integrated stress response. ATF6, also metazoan-specific and with 
ot [GenBank:NM_007348.3] and p 1 [GenBank:NM_017406.4] 
paralogs, is resident to the ER but transits to the Golgi during 
stress, where it is cleaved by regulated intramembrane proteolysis 
to liberate an active bZIP transcriptional activator. Together, these 
bZIPs coordinate enhancement of protein synthesis, degradation, 
folding, modification, and trafficking through gene regulation. 

An oft overlooked feature of UPR activation is that between 
20 and 50 percent of regulated genes are actually suppressed 
by ER stress depending on the conditions, yet much less is 
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known about the mechanisms responsible for this suppression 
and the physiological consequences thereof. A portion of this 
suppression can be attributed to regulated IREl-dependent 
decay, in which the IRE1 endonuclease degrades ER-associated 
mRNAs (Hollien and Weissman, 2006). Transcriptional mech- 
anisms for suppression have been identified as well, including 
direct suppression by the bZIP C/EBP family member CHOP 
[GenBank:NM_007837.3] (Ron and Habener, 1992), titration of 
the coactivator CRTC2 [GenBank:NM_028881.2] (Wang et al, 
2009), and translational regulation of the suppressive LIP isoform 
of C/EBPF5 [GenBank:NM_009883.3] (Li et al, 2008; Arensdorf 
and Rutkowski, 2013). Each of these mechanisms was identified 
through the behavior of target genes, so the extent to which any 
of them contributes to global gene suppression is not clear. 

In the liver, the most evident consequence of ER stress is 
lipid accumulation (Rutkowski et al., 2008; Yamamoto et al., 
2010; Zhang et al, 2011). This lipid accumulation, or steatosis, is 
accompanied by suppression of a host of genes involved in hepatic 
lipid metabolic processes, including fatty acid oxidation, lipo- 
genesis, cholesterologenesis, and VLDL production. Given that 
some of these processes are mutually antagonistic (e.g., fatty acid 
oxidation and lipogenesis), it seems likely that some are sup- 
pressed as primary responses to ER stress, and others as secondary 
consequences of feedback mechanisms. 

Some regulation of hepatic lipid metabolism can be attributed 
directly to the action of canonical UPR signaling. XBP1 
can bind to the promoters and stimulate transcription of 
lipogenic genes (Lee et al., 2008) and of the ER oxidore- 
ductase PDI [GenBank:NM_001032.2] (Wang et al, 2012), 
the latter of which stimulates VLDL secretion by virtue of 
its interaction with the microsomal triglyceride transfer pro- 
tein (MTTP) [GenBank:NM_001163457.1]. IREla can also 
directly degrade mRNAs encoding lipogenic genes when Xbpl 
is ablated through its regulated IREl-dependent decay activ- 
ity (So et al, 2012), thus acting at cross-purposes with its 
downstream target XBPL The ER-localized transcription fac- 
tor CREBH [GenBank:NM_145365.3] also contributes to lipid 
homeostasis (Zhang et al, 2012), although it is not yet clear 
whether that function is direct or indirect. 

Steatosis arises as a common phenotype in response to ER 
stress when any of the UPR signaling pathways is ablated, or when 
UPR signaling is intact but ER protein folding is compromised by 
deletion of the ER cochaperone p58 IPK [GenBank:NM_008929.3] 
(Rutkowski et al., 2008). This steatosis is likely at least partially 
caused by impaired secretion of triglyceride-rich VLDL parti- 
cles from the stressed ER (Ota et al., 2008; Rutkowski et al., 
2008; Caviglia et al., 2011) and enhanced uptake (Jo et al., 2013). 
However, ER stress also elicits extensive alterations in the expres- 
sion of genes involved in lipid metabolism, and these alterations 
are more severe and persistent when any branch of the UPR — 
or when ER protein folding — is compromised (Rutkowski et al., 
2008). Thus, these alterations correlate with the development of 
steatosis, although it is not known which events precede lipid 
accumulation and which follow as a consequence. That they 
emerge irrespective of which UPR pathway is ablated argues 
that most such metabolic genes are not directly regulated by 
ATF4, ATF6, or XBP1, but by some mechanism or mechanisms 



that indirectly tie metabolic gene regulation to the ER stress 
burden. To some degree the dysregulation of lipid metabolism 
can be attributed to CHOP (Chikka et al., 2013), which is 
redundantly regulated by each of the three UPR pathways and 
which is expressed more robustly when stress is more severe — 
as when the UPR or ER protein folding is disrupted (Rutkowski 
et al., 2008). However, Chop~/~ animals are only partially pro- 
tected from hepatic steatosis during ER stress (Rutkowski et al, 
2008), suggesting that other as yet uncovered pathways exist 
as well. 

Given the extensive nature of metabolic gene regulation during 
ER stress, there likely exists a mechanistic hierarchy of regulation, 
with some metabolic genes being more proximally connected 
to UPR pathways and others lying downstream of these initial 
events. However, the global organization of lipid metabolic gene 
regulation during ER stress has not been studied. Thus, our goal 
in this work was to begin to decipher the structure of ER stress- 
mediated metabolic gene regulation by establishing the temporal 
progression of such events in the mouse liver, and to infer hierar- 
chical relationships using a functional genomics approach, based 
upon the behavior of coordinately regulated groups of genes in 
wild-type mice vs. mice lacking the ER stress sensor ATF6a. 

MATERIALS AND METHODS 
ANIMAL EXPERIMENTS 

All protocols for animal use were reviewed and approved by 
the University Committee on Use and Care of Animals at the 
University of Iowa or the University of Michigan. Animals were 
bred in house, and were fed standard rodent chow and housed 
in a controlled pathogen free environment with 12 h light and 
dark cycles. Animals used were of varying ages and genders, with 
control and experimental groups having similar composition. 
Animals were fasted for 4 h prior to sacrifice, which was carried 
out in daytime hours. 

LIPID ANALYSIS 

ADRP immunostaining was as described (Chikka et al., 2013). 
For the trigylceride assay, a 100 mg piece of liver was homog- 
enized in 1 ml of ice cold extraction buffer ( 1 mM Tris pH 
7.6, 1 mM EGTA, 1 mM MgC^) containing protease inhibitors. 
A 200 [iL aliquot of the homogenate was placed into a new 
1.5 ml tube on ice, and an additional aliquot was set aside to 
determine the protein concentration of each sample. A 750 \iL 
aliquot of a Chloroform:Methanol mixture (1:2 ratio) was added 
to the 200 [lL homogenate sample and vortexed vigorously for 
15 s. The samples were incubated at room temperature for 
1 h, with the samples vortexed every 15min. Following the 
hour incubation, 250 |iL of Chloroform was added to each 
sample and vortexed for 15 s, then incubated at room tem- 
perature for 15min. Two hundred u,L of distilled water was 
then added to the sample and vortexed as above. The sam- 
ples were centrifuged at 5000 rpm for lOmin, and the bottom 
organic layer was collected and placed in a fresh tube. The 
sample was evaporated under nitrogen gas, and the remain- 
ing lipids were dissolved in 200 [iL of isopropanol. From 
these samples, the triglyceride levels were determined using the 
Infinity Triglycerides Reagent (Thermo Scientific, TR22421) per 
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manufacturer's instructions. The Cayman triglyceride standard 
(Cayman Chemical, 10010509) was used to generate a standard 
curve. Oil Red O staining was as described (Rutkowski et al., 
2008). 

RNA ANALYSIS 

The 8 h microarray has been published (Rutkowski et al, 2008). 
For the 34 h microarray, mice were injected with 1 mg/kg TM, 
and mRNA was prepared from isolated livers and analyzed by 
Affymetrix microarray in exactly the same way. The NCBI GEO 
accession number for both arrays is GSE48939. Expression cat- 
egories for each probeset are provided in Table SI. RT-PCR and 
qRT-PCR analysis, including validation of all primer sets, was as 
previously described (Rutkowski et al, 2006; Tyra et al., 2012), 
except that gene expression was normalized against the average 
expression of two housekeeping genes (Btf3 and Ppia) rather 
than only one. Primer sequences can be found in (Rutkowski 
et al., 2006, 2008; Wu et al, 2007; Tyra et al, 2012) and 
Table S2. 

BI0INF0RMATIC ANALYSIS 

Pathway enrichment analysis was performed using FunNet soft- 
ware (Prifti et al., 2008). Data represent GO biological process 
annotations with a decorrelated enrichment computation and a 
5% false discovery rate correction. Transcription factor binding 
site analysis was performed using oPOSSUM software (Ho Sui 
et al., 2005). Data represent a single site analysis of vertebrate 
transcription factor binding sites within 2000 base pairs upstream 
or downstream of the transcription start site. The results of 
this analysis were visualized using Cytoscape software. The net- 
work was limited to genes with GO annotations involved in lipid 
metabolism and transport using the BiNGO plug- in (Maere et al., 
2005). 

CHROMATIN IMMUN0PRECIPITATI0N 

Livers were isolated from 6-8 week-old mice and prepared for 
ChIP using the ChIP Tissue Chromatin Shearing Kit with SDS 
(Covaris). Samples were sonicated using a S220 focused ultrason- 
icator (Covaris) to produce DNA bands between 100 and 1000 bp. 
Following sonication, the immunoprecipiation was carried out 
as described in (Arensdorf and Rutkowski, 2013) using HNF4a 
antiserum (H-171, Santa Cruz) or non-specific IgG (12-370, 
Millipore). 

RESULTS AND DISCUSSION 

SUPPRESSION OF METABOLIC TRANSCRIPTIONAL MASTER 
REGULATORS PRECEDES LIPID ACCUMULATION DURING ER STRESS 

Lipid accumulation can be induced in the liver by exposure 
to the inhibitor of N-linked glycosylation tunicamycin (TM) 
or the proteasomal inhibitor bortezomib, or by overexpression 
of a misfolded ER client protein such as coagulation Factor 
VIII [GenBank:NM_001161373.1] (Rutkowski et al, 2008; Zhang 
et al., 2011; Chikka et al., 2013). Each of these treatments acti- 
vates either the UPR or integrated stress response and each leads 
to qualitiatively similar changes in the expression of key metabolic 
genes in the liver. These genes include both transcription factors 
and cofactors involved in controlling metabolism as well as the 



downstream targets of these factors that encode the key func- 
tional enzymes involved in lipid catabolism, anabolism, storage, 
and secretion (Rutkowski et al, 2008). 

The temporal organization of these gene regulatory changes 
prior and subsequent to the onset of lipid accumulation has not 
been analyzed. It is likely that the earliest of these regulated events 
are directly mechanistically connected to the UPR and/or inte- 
grated stress response. Thus, we examined the time course of 
hepatic lipid accumulation in wild-type mice in response to TM, 
which is a more robust inducer of ER stress and steatosis than 
other stimuli (Chikka et al., 2013). We monitored lipid accumu- 
lation by three distinct criteria: accumulation of Oil Red O in 
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FIGURE 1 | ER stress causes substantial hepatic lipid accumulation 
within 8 h. (A) C57BL/6J mice were challenged with 1 mg/kg TM for the 
indicated times, livers were frozen in OCX and lipids were stained with Oil 
Red O. Scale bar = 50 |im. (B) Same as (A), except lipid content was 
assessed by immunohistochemical staining for the lipid droplet marker 
protein ADRR Scale bar = 50 |im. (C) Same as (A), except triglyceride 
content was measured by colorimetric assay after extraction of neutral 
lipids, p < 0.001 by One-Way ANOVA. n = 3 samples per time point. Error 
bars here and elsewhere show means ± SDM. 
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lipid droplets of fresh frozen liver sections (Figure 1A); immuno- 
histochemical detection of the lipid-droplet associated pro- 
tein adipophilin (ADRP) [GenBank:NM_007408.3] (Figure IB); 
and direct biochemical assessment of hepatic triglyceride levels 
(Figure 1C). Results from all three assays were similar: hepatic 
lipid content increased most substantially at the 8 h time point. 
Therefore, the key genetic regulatory events responsible for alter- 
ing lipid metabolism are likely to occur prior to 8 h, while those 
that occur later than 8 h after challenge are more likely secondary 
effects that contribute to lipid disruption tangentially, if at all. 

Next, we examined the timing of UPR activation and of 
the regulation of metabolic genes. TM led to maximal IRE la- 
dependent splicing of Xbpl mRNA within 2 h; this splicing 
persisted through 8h but was diminished at later time points 
(Figure 2A). Likewise, every UPR- regulated target gene was 
significantly upregulated by 2 h, peaked at 4-8 h, and diminished 



thereafter (Figure 2B). These genes depend to varying extents 
on activity of each of the three UPR pathways (Harding et al., 
2003; Lee et al, 2003; Wu et al, 2007). That they are uni- 
formly upregulated by 2 h suggests that each of the three canoni- 
cal UPR-regulated transcriptional activators — ATF6a, ATF4, and 
XBP1 — is functionally active by this very early time point. 

We then analyzed expression of genes involved in lipid 
metabolism, including both transcriptional regulators and key 
rate-limiting metabolic enzymes. We grouped the genes accord- 
ing to the time point at which they were first down-regulated. 
We immediately observed that, in contrast to conventional UPR 
target genes, the regulation of metabolic genes occurred in stages. 
The earliest group regulated included select transcription factors 
and cofactors. The coregulator Pgclb [GenBank:NM_133249.2] 
was downregulated by 2 h although not at 4 h (Figure 2C), mak- 
ing the true timing of its suppression ambiguous. In contrast, by 
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FIGURE 2 | Staggered suppression of metabolic genes during ER stress. 

(A) The spliced (spl) and unspliced (us) forms of Xbpl mRNA were detected 
by RT-PCR of total RNA isolated from the livers of mice treated with vehicle 
(veh) or 1 mg/kg TM for the indicated times. Each lane shows a separate 
animal. Image is shown in black-to-white inverted form for greater visual 
clarity. (B) Expression of the indicated UPR target genes was determined by 
qRT-PCR from the animals shown in (A), with Btf3 and Ppia used as 



normalizing controls. Expression here and in subsequent figures is given on a 
log2 scale relative to the vehicle-treated condition. Here and elsewhere 
unless noted: *p < 0.05; *p < 0.1 by two-tailed student's f-test. (C-F) 
Expression of metabolic genes was assessed by qRT-PCR as in (B), and 
genes were grouped according to the time point at which downregulation 
(p < 0.1 ) was first observed. The process in which each gene participates is 
listed. 
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4 h the coregulator Pgcla [GenBank:NM_008904.2] and the tran- 
scriptional activators Cebpa [GenBank:NM_007678.3] and Ppara 
[GenBank:NM_001 113418.1] were suppressed (Figure 2C). 
PGC-ip contributes to lipogenesis, fatty acid oxidation, and 



A vehicle 
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VLDL production and secretion (Lee et al., 2003; Lin et al., 2003, 
2005a,b; Wolfrum and Stoffel, 2006), and PGC-la contributes 
to the latter two of these processes (Louet et al., 2002; Lin et al., 
2005a; Rhee et al., 2006). C/EBPa has general roles in energy 
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FIGURE 3 | Loss of Atf6a exacerbates steatosis and long-term metabolic 
gene suppression. (A) Wild-type or Atf6a.~l~ mice were challenged with 
1 mg/kg TM or vehicle for 48 h, and ADRP immunostaining was carried out as 
in Figure 1. Scale bar= 50|j,m. (B-F) Wild-type or Atf6a~t~ mice were 
challenged with 1 mg/kg TM or vehicle for 34 h, and global mRNA expression 
was assessed by Affymetrix microarray. Array-determined expression of the 



indicated genes is shown. Genes were arranged in the same groupings as in 
Figure 2. Statistical significance was calculated by two-tailed student's f-test, 
comparing expression in TM-treated Atf6arl~ animals against TM-treated 
wild-type animals, n = 3 animals per group. (G) Mice were treated with TM 
or vehicle for 48 h as in (A), and expression of the indicated genes was 
determined by qRT-PCR. Significance was determined as in (B-F). 
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homeostasis including lipid and glucose metabolism, while 
PPARct is a master regulator of fatty acid oxidation (Desvergne 
et al, 2006). 

The genes downregulated at later times — i.e., after the onset 
of pronounced lipid accumulation — included both transcrip- 
tional regulators and downstream genes encoding rate-limiting 
enzymes in metabolic processes, as indicated in Figures 2D-F. 
These results suggest that PGC-lct, C/EBPot, and PPARa (and 
possibly PGC-1{5) are most likely the first lipid metabolic genes 
regulated by ER stress. Further, the processes downstream of these 
factors — namely, fatty acid oxidation and lipoprotein biogenesis 
and transport — are the first affected via gene regulatory mecha- 
nisms. In addition, because lipogenic genes are not altered until 
much later (18 h), they suggest that inhibition of lipogenesis is a 
secondary consequence of ER stress, perhaps occurring as a result 
of negative feedback when lipids begin to accumulate. In addi- 
tion, even though XBP1 splicing is induced by ER stress, we found 
no evidence that lipogenic gene expression was stimulated under 
these conditions. 

DELETION OF ATF6« EXACERBATES METABOLIC GENE SUPPRESSION 

Mice lacking ATF6ct, while otherwise apparently normal and 
healthy, exhibit dramatic steatosis upon challenge with TM 
(Rutkowski et al, 2008; Yamamoto et al., 2010). This phenotype is 
illustrated by ADRP staining in Figure 3A; it arises from an inabil- 
ity to restore ER homeostasis upon challenge. Therefore, we rea- 
soned that these mice could be used to expose the stress-regulated 
gene expression changes that truly underlie lipid dysregulation, 
since those changes should be amplified in Atf6a~l~ animals. A 
second benefit of these animals is that they can be used to identify 
truly stress-responsive expression changes; those that result from 
ER stress-independent properties of TM would not be expected 
to differ between wild-type and Atf6arl~ animals, since ATF6a 
is not thought to act outside the context of ER stress, and its dele- 
tion does not alter the apparent pharmacological activity of TM 
(Rutkowski et al, 2008). 

We approached this goal by challenging wild-type and 
Atf6aT^ animals for an extended period (34 h) with TM, 
and then profiling global gene expression by microarray. We 
chose this approach in part because we had previously used 
microarray profiling to characterize gene expression in these 
same mouse strains following 8 h of TM challenge (Rutkowski 
et al, 2008). Conducting a second microarray study with an 
identical gene chip at a later time point would provide us 
with the unique opportunity to determine, for every gene 
on the array, its expression in two different genotypes, under 
two different treatment regimens (vehicle and TM) at two 
different times. 

By having 8 distinct combinations of experimental condi- 
tions, we anticipated that we could begin to identify clusters of 
coordinately regulated genes, including those metabolic genes 
most responsive to ER stress. This relied upon the assump- 
tion that genes which are part of a common functional path- 
way (in this case, lipid metabolism) should be regulated simi- 
larly. This approach has been used to great effect in yeast, and 
more recently in cultured mammalian cells, to expose previ- 
ously hidden components of the secretory apparatus and other 



pathways of interest (Schuldiner and Weissman, 2013). Here, our 
goal was to infer hidden transcriptional regulators using temporal 
patterns of gene expression as an output. 

From the 34 h microarray, we examined the expres- 
sion of the UPR and metabolic genes shown in Figure 2. 
Consistent with the role of ATF6ct in contributing to chap- 
erone expression (Wu et al, 2007; Yamamoto et al., 2007), 
most ER chaperones and cochaperones known to be regulated 
by ATF6a, such as DerB [GenBank:NM_024440.2], Erdj3 
[GenBank:NM_001 190804.1], Grp94 [GenBank:NM_0 11631.1], 
p58 IPK , and Erp72 [GenBank:NM_009787.2], were not upregu- 
lated to as great an extent in Af/6a~/~ animals as in wild-type 
animals (Figure 3B). Conversely, expression of Gadd34 
[GenBank:NM_008654.2], which largely depends upon the 
PERK axis of the UPR (Marciniak et al, 2004), was elevated 
in Aif6oi~/~ mice, consistent with persistent ER stress and 
activation of the other limbs of the UPR. 

Amongst the metabolic genes examined, none was 
upregulated by ER stress in either genotype at 34 h. 
Two genes — the lipogenic and cholesterologenic reg- 
ulators Srebfl [GenBank:NM_01 1480.3] and Srebf2 
[GenBank:NM_033218.1] — were equally downregulated in 
both genotypes (Figures 3D,E). However, the large majority 
of genes were expressed at normal or near-normal levels in 
wild-type animals, but deeply suppressed in Af/6a~/~ animals 
(Figures 3C-F). qRT-PCR profiling of a sampling of these genes 
from an independent experiment confirmed these findings 
(Figure 3G). These data are consistent with the idea that acute 
ER stress inhibits fatty acid oxidation and lipoprotein biogenesis 
at the level of gene expression, and does not stimulate lipogenic 
gene expression. 

FUNCTIONALLY RELATED GENE GROUPS CLUSTER TEMPORALLY 

Having these 8 distinct experimental conditions enabled us 
to compare the global behavior of hepatic gene expression in 
response to ER stress at early vs. late times in normal ani- 
mals vs. those with a compromised UPR. A heatmap showing 
TM-regulated genes makes two observations clear: First, gene 
expression differences between the two genotypes were far more 
extensive at 34 h than at 8 h. Second, at 34 h, many genes have 
returned to normal or near-normal expression in wild-type ani- 
mals, but have remained regulated, or become even more regu- 
lated in the same direction, in At/6a~/~ animals. 

Each gene was then categorized based on whether it was upreg- 
ulated, downregulated, or unchanged by ER stress in wild-type 
and Af/6a~/~ animals after 8 or 34 h of stress. Thus, a gene 
could fall into any of 81 (3 x 3 x 3 x 3) expression profiles. The 
assignments for all probesets are provided in Table SI. We were 
able to analyze the genes in this way in part because very few 
genes differed in their basal (i.e., unstressed) expression between 
genotypes; almost all genotype-dependent changes in expression 
were caused by TM, so basal expression differences were not 
confounding (Rutkowski et al, 2008). 

Because Atf6aT^ animals became more steatotic than wild- 
type animals, genes that were not differentially expressed between 
the two genotypes upon 34 h of TM treatment were not pursued 
further, as these were unlikely to be major contributors to the 
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FIGURE 4 | Clustering of genes according to temporal regulation in 
wild-type and Atf6a~l~ animals. (A) The expression of every probeset on the 
Affymetrix microarray described in Figure 3 was aggregated with expression 
data from a previously published identical array comparing gene expression in 
wild-type and AtfQarl~ animals 8 h after challenge with vehicle or 2 mg/kg TM 
(Rutkowski et al., 2008). For each time point, expression was determined using 
a log2 scale relative to vehicle-treated wild-type animals at that time point. Only 
probesets showing significant (p < 0.05) expression differences (1.5-fold, or 
±0.58 on the log2 scale) in one or more of the four of the four conditions 
(wild-type or/4ff6rW~ at 8 or 34 h) are shown by heatmap, which accounted for 
~7000 of the ~45,000 probesets on the array. Extent of up- or downregulation 
is shown by intensity of red or blue coloration, respectively. Each column 
depicts expression level averaged among the three animals per group. (B,C) 
Every probeset on the array was characterized by its expression in the following 
four ways, with differences defined as > 1.5-fold, p < 0.05: (1 ) up, down, or 



unchanged in wild-type TM-treated animals at 8 h relative to vehicle-treated 
wild-type; (2) up, down, or unchanged in Atf6a~t~ TM-treated animals at 8 h 
relative to TM-treated wild-type; (3) same as (1) but 34 h; and (4) same as (2) but 
34 h. The genes that showed a difference by criterion (4) were broken down into 
groups based on their behavior with respect to these criteria, and the number of 
genes in the nine most populated groups is shown in (B). The group of genes 
that were upregulated by ER stress in wild-type animals at both time points, but 
were less upregulated in AtfQarl~ animals — i.e., genes that could be 
understood as directly ATF6a-dependent — is also accounted for. (C) provides a 
key for illustration of gene expression patterns. (D-G) Expression pattern for 
each of the gene groups shown in (B). These include genes shown in Figure 2 
(those that did not fall into one of these groups are illustrated in Figure SI) as 
well as genes involved in ER protein processing found in Group E. For genes 
represented by more than one probeset, the behavior most commonly 
represented and/or most consistent with qRT-PCR data is shown. 



steatotic phenotype. Although 54 remaining combinations of 
expression were possible (3x3 x 3x2), 90 percent of the 
genes fell into only 9 categories (A-I, Figure 4B). Each cate- 
gory of genes was then given a graphical representation of the 
expression pattern of its members (Figures 4C-G). The two 
most populated groups of genes (Groups A and B) were those 
that were unaffected by stress at 8 h in either genotype, or in 
wild-type animals at 34 h, but were up- or down-regulated, 
respectively, at 34h in Atf6a~'~ animals (Figure4E). Among 
the metabolic genes depicted in Figures 2, 3, a number of 
these populate Group B, and they include genes encoding rate- 
limiting enzymes in each pathway of fatty acid oxidation — Acoxl, 
Cptla, and Cyp4al0 [GenBank:NM_015729.3, NM_013495.2, 
and NM_0 100 11.3], which control peroxisomal, mitochondrial, 
and microsomal oxidation, respectively. As the expression of these 
genes is not altered until the later timepoint, they are unlikely to 
be proximally connected to the UPR, but more likely represent 



indirect effects — for example, of suppression of PPARa. The 
next two most populated groups (Groups C and D) included 
those genes that were either up- or down-regulated early in 
both genotypes, and whose expression returned to normal lev- 
els in wild-type animals by 34 h but which remained regulated in 
Atf 6a~/~ animals. Among metabolic genes, Group D included 
the coregulators Pgcla and Pgclb. Closely related to these groups 
were Groups E and F, which included genes that were up- or 
downregulated early in both genotypes, and for which this reg- 
ulation was enhanced in Atf6aT^ animals at 34 h. The other two 
rapidly regulated metabolic genes, Cebpa and Ppara, were found 
in Group F. 

The lipogenic genes Fasn [GenBank:NM_007988.3] and Acaca 
[GenBank:NM_133360.2] and the cholesterologenic genes Acat2 
[GenBank:NM_009338.3] and Hmgcr [GenBank:NM_008255.2] 
showed no evidence of upregulation in the 8 or 34 h array data, 
and Fasn and Acat2 were actually suppressed, as were Srebfl and 
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Srebf2 themselves (Figure SI). While lipogenesis can be stim- 
ulated by non-transcriptional mechanisms — most notably pro- 
cessing of SREBP-lc and SREBP-2 — this processing would result 
in stimulation of their downstream target genes, which is not 
evidenced here. 

Groups D and F were of the most interest to us because 
they represented those genes most likely to be proximally mech- 
anistically connected to UPR signaling, since they were regu- 
lated rapidly by ER stress and since their long-term expression 
coincided with the persistent ER stress and worsening lipid 
accumulation seen in Atf6a~/~ animals. Group F in particular 
was noteworthy because its counterpart cohort of upregulated 
genes — Group E — included a number of genes known to be 
direct targets of UPR transcription factors and which are thus 
proximally connected to UPR signaling (Figure 2D). 

Also supporting the validity of the approach was the group 
of genes that were upregulated by ER stress in wild-type 
animals at both time points, but that were not as upregu- 
lated in Atf6a~^~ animals at both time points (Figure 4G). 
These genes would fit the expected profile of direct targets of 
ATF6a. While there were relatively few genes in this group, they 
included those already described as direct ATF6a targets, includ- 
ing Erp72, p58 IPK , Erdj3, and Derl3 (Wu et al, 2007; Yamamoto 
et al., 2007). This finding supports the idea that coregulated 
genes can be discriminated based on their expression profile 
in this setup. 

This idea was reinforced by pathway analysis. The genes 
from each of the ten groups (A-I and ATF6) were analyzed for 
functional enrichments of Gene Ontology (GO) pathways using 
FunNet (Prifti et al., 2008). As proof-of-concept, pathway analy- 
sis of Group ATF6 genes yielded "unfolded protein response" as 
the most significantly enriched process, which would be expected 
of a group encompassing ATF6a direct targets (Figure 5A). Genes 
representing processes relevant to lipid metabolism were enriched 
in all of the downregulated groups — B, D, F, and G — but not 
in the upregulated groups A, C, E, H, and I (Figures 5B,C). 
Conversely, each of the upregulated groups was enriched in genes 
representing pathways relevant to protein synthesis, trafficking, 
and degradation. UPR activation during ER stress is known to 
transcriptionally augment the cellular protein biogenesis machin- 
ery through the action of the major UPR-regulated transcrip- 
tional activators XBP1, ATF4, and ATF6 (Harding et al, 2003; Lee 
et al., 2003; Wu et al., 2007). This pathway analysis thus suggests 
that, at least in the liver, suppression of genes involved in lipid 
metabolism represents a concerted focus of UPR activation. 

FUNCTIONAL GENOMIC ANALYSIS IDENTIFIES HNF4a AS A PROXIMAL 
REGULATOR OF METABOLIC GENE EXPRESSION DURING ER STRESS 

We wished to harness the statistical power of our microar- 
ray comparisons in order to predict regulatory transcription 
factors whose activity was altered by ER stress; these would 
be the most likely to be proximally mechanistically connected 
to the UPR. To accomplish this, we subjected the genes in 
each group to analysis using oPOSSUM (Ho Sui et al., 2005), 
which searches the promoter/enhancer region of each gene for 
consensus transcription factor binding sites from the JASPAR 
CORE database. Establishing the validity of this approach, the 



genes in Group ATF6 yielded the transcription factor NFYA 
[GenBank:NM_001 110832.1] as the sole statistically significant 
hit (Figure 6A). ATF6a is known to dimerize with NFYA to regu- 
late transcription from promoters containing ERSE and ERSE II 
elements (Yoshida et al, 2000, 2001; Kokame et al, 2001). NFYA 
was not enriched in any other group, underscoring the specificity 
of the algorithm. 

Because genes involved in lipid metabolism were found in 
groups B, D, F, and G, we sought hits found in those groups and 
no others, with a particular emphasis on groups D and F, since 
these groups contain the early-responding genes. Remarkably, 
binding sites for two transcription factors — HNF4a and NR2F1 
[GenBank:NM_010151.2] — were enriched in Groups B, D, and 
F but not in any other group (Figures 6B,C). These factors 
have similar consensus binding sites (Kimura et al., 1993; Ellrott 
et al., 2002; Schmidt et al., 2010), explaining why they segregate 
together in this analysis. This finding suggests that the activ- 
ity of these two transcriptional activators is altered during ER 
stress, since genes containing binding sites for these factors are 
downregulated by ER stress, and their downregulation is exacer- 
bated by the ongoing stress seen in Atf 6a~/~ mice. Conversely, 
groups A, C, and E were enriched for genes with potential bind- 
ing sites for ELK4, also known as SRF Accessory Protein (SAP)-l 
[GenBank:NM_007923.2]. 

Group I, encompassing genes upregulated by long-term 
ER stress in wild-type animals and further upregulated in 
Atf6oT^ animals, was enriched for binding by CREB1 
[GenBank:NM_001037726.1], which shares a consensus bind- 
ing sequence with the ER-localized transcription factor CREBH 
(Zhang et al, 2006). CREBH is activated by proteolysis induced 
by, among other stimuli, ER stress, and it regulates expression 
of acute phase response genes and metabolic genes (Zhang et al, 
2006, 2012; Luebke -Wheeler et al, 2008). Accordingly, the genes 
in Group I have a profile expected for CREBH targets — namely 
that they are upregulated by ER stress, and further upregulated in 
Af/6a~/~ animals, in which ER stress is exacerbated. However, 
lipid metabolic genes were not enriched in Group I, and the 
CREB1 binding site was not enriched in the groups containing 
lipid metabolic genes. In addition, other than the gene encoding 
serum amyloid P-component (Apes) [GenBank:NM_011318.2], 
most genes encoding acute phase response proteins (SAA pro- 
teins, CRP, coagulation and clotting proteins, etc.) were not 
significantly upregulated at either time point in either geno- 
type, suggesting that CREBH activity during bona fide ER stress 
(as opposed to endotoxin, pro-inflammatory cytokines, or other 
stimuli) might be minimal, and the genes in Group I might 
instead be regulated by another factor from the CREB family. 

ELK4 is part of the Ternary Complex Factor (TCF) fam- 
ily of transcription factors that interact with Serum Response 
Factor (SRF) (Dalton and Treisman, 1992). A role for ELK4 in 
hepatic gene expression has not been described, nor has ELK4 
been directly linked to ER stress. ELK4 has been shown to be 
activated by JNK-dependent phosphorylation (Janknecht and 
Hunter, 1997), and JNK [GenBank:NM_016700.4] is activated 
by ER stress (Urano et al., 2000). Thus, we speculate that ELK4 
activity might be regulated during ER stress in the liver by JNK or 
other MAP kinases to promote expression of genes in groups A, 
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FIGURE 5 | Functionally related genes cluster by temporal regulation. 
(A-C) Each of the gene groups in Figure 4 was subjected to Gene Ontology 
pathway analysis using FunNet. The top seven most significant pathway 
enrichments were then reordered by the number of genes from that group 
that were pathway "hits," with pathways having the most "hits" listed higher. 



C, and E. NR2F1, also known as COUP-TF, is an orphan receptor. 
Deletion leads to perinatal lethality with extensive dysregulation 
of neuronal differentiation (Qiu et al., 1997). Its function in 
the liver is less clear, although it has been shown to coactivate 
transcription synergistically with HNF4a (Ktistaki and Talianidis, 
1997; Yanai et al., 1999). 

HNF4ct is expressed most strongly in the liver, intestine, kid- 
ney, and pancreas. Deletion is lethal during gastrulation (Chen 
et al., 1994). Mice with a liver-specific deletion of HNF4a develop 
steatosis concomitant with impaired ApoB and Mttp expres- 
sion and VLDL production (Hayhurst et al., 2001). Hepatic 
knockdown of HNF4a by adenoviral delivery resulted in steato- 
sis and in impaired VLDL production, along with essentially 
uniformly diminished expression of a host of genes involved 
in lipid metabolism, including many of those reported here 
(Yin et al., 2011). Thus, loss of HNF4a phenocopies many of 
the lipid metabolic genetic changes that are induced by ER 
stress. Overexpression of HNF4ot in primary hepatocytes resulted 
in upregulation of many of these same genes, though not of 
Srebfl nor Srebf2 (Yin et al, 2011). This exception is notable 
because Srebfl and Srebf2 are conspicuous among the metabolic 
genes analyzed here in the fact that their downregulation is not 
exacerbated at 34h in Atf6aT^ mice (Figures 3D,E). HNF4a 



For reasons of space, five of these seven pathways from each group are 
shown. In no case was a lipid metabolism process enriched among the 
upreguiated gene groups. In (B), pathways relevant to lipid metabolism are 
listed in black, and other pathways in gray. In (C), pathways relevant to 
protein processing are listed in black. 



binding sites in the mouse liver genome have been characterized 
by ChlP-seq (Schmidt et al., 2010), and genes in Groups B, D, and 
F are confirmed HNF4a targets (Figure 6D). 

Our results, together with the known activity of HNF4a, pre- 
dict that ER stress leads to diminished activity of HNF4a, and that 
this occurs as an early event in metabolic gene regulation by ER 
stress. To test this prediction, we analyzed the binding of HNF4a 
to ChlP-seq-defined sites in the promoter/enhancer regions of 
three of the four earliest regulated metabolic genes — the tran- 
scription regulators Cebpa, Pgcla, and Ppara (the proximal Pgclb 
promoter/enhancer did not contain a predicted or validated 
HNF4a binding site). We found that treatment of animals with 
TM for 8 h did not change either the total expression (Figure 7A) 
or nuclear localization (Figure 7B) of HNF4a. Yet, consistent 
with our prediction, we found that HNF4a chromatin binding 
decreased significantly at several sites in the Cebpa and Pgcla 
promoters upon treatment of animals with TM (Figure 7C). 
This diminishment was not uniform; several sites within the 
three promoters showed unaltered HNF4a binding (Figure 7A). 
Together, these findings suggest that ER stress reduces the 
activity of HNF4a at specific sites in the genome through 
a mechanism that is independent of the HNF4ct expression 
level. 
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FIGURE 6 | Transcription factor prediction implicates ELK4, HNFIce, 
NR2F1, and HNF4« as hidden regulatory nodes in hepatic 
stress-dependent gene regulation. (A-C) Each of the gene groups in 
Figure 4 was subjected to oPOSSUM single-site analysis, which searches 
regulatory regions (in this case, ±2000 bp from the transcriptional start site) 
for potential binding sites of transcription factors identified in the JASPAR 



CORE database. The results were limited to a Z-score > 10 and a Fisher 
score < 0.01. (D) The data from (C) were visualized using the BINGO plug-in 
for Cytoscape software, considering only genes relevant to lipid metabolism 
as annotated from GO analysis. oPOSSUM-predicted binding sites are shown 
using dashed lines, while genes with confirmed HNF4a-binding sites are 
shown by solid lines. 



Frontiers in Genetics | Genomic Endocrinology 



September 2013 | Volume 4 | Article 188 | 10 



Arensdorf et al. 



Gene suppression by the UPR 



CONCLUSIONS 

The hepatic lipid dysregulation that is elicited by ER stress is 
accompanied by sweeping alterations to the expression of genes 
involved in the process. With few exceptions, these genes are 
downregulated by ER stress, and encompass numerous metabolic 
pathways including fatty acid oxidation, lipogenesis, triglyceride 
storage and secretion, and phospholipid synthesis. The down- 
regulated genes include both those encoding the key enzymes of 
each of these processes as well as the upstream transcription fac- 
tors that regulate them. However, it would be unlikely that these 
genes would be directly suppressed by canonical UPR-regulated 
transcription factors, since these are, for the most part, tran- 
scriptional activators rather than repressors. Further, ablation of 
the ATF6ol pathway of the UPR compromises recovery from ER 
stress and leads to an exacerbated steatotic phenotype, yet ATF6a 
does not directly act upon genes involved in lipid metabolism — 
no such genes populate either group ATF6 or its downregulated 
converse. Therefore, other mechanisms linking UPR activation to 
metabolic gene regulation must be at work. The extensive nature 
of metabolic gene regulation during ER stress also suggests a hier- 
archical organization of gene regulatory events, with most genes 
regulated as an indirect consequence of earlier events. Such an 
organization predicted that the temporal ordering of gene regu- 
lation could be used to identify the transcriptional events most 
proximal to UPR signaling. 
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FIGURE 7 I Diminished HNF4a binding at Cebpa and Pgda 
promoters during ER stress. (A,B) Wild-type mice were challenged 
with 1 mg/kg TM for 8h, and expression of HNF4a was assessed by 
immunoblot (A) or immunohistochemistry (B). Scale bar = 50|j,m. (C) 
HNF4a binding to the regulatory regions of the indicated genes was 



With these facts in mind, we hypothesized that a "bottom-up" 
approach could be used to establish which ER stress-dependent 
metabolic gene expression changes occurred earliest, and to 
identify common regulators of those genes. To that end, our 
analysis predicted altered activity of HNF4a, which we then 
confirmed experimentally. We can conclude that its diminished 
binding to the promoters of Cebpa and Pgcla is likely to con- 
tribute to the suppression of these genes, because knockdown 
of HNF4a is already known to have similar effects on metabolic 
genes to those reported here (Yin et al., 201 1). We have combined 
the experimental and bioinformatic experiments described here 
with existing literature on the roles of HNF4a, C/EBPa, PGC-la, 
and PPARot into a working model describing the genetic hierar- 
chy of lipid metabolic gene regulation during stress (Figure 8). 
Although these relationships were elicited using TM, the observa- 
tion that proteasome inhibition or overexpression of a misfolded 
secretory protein leads to similar genetic changes and lipid accu- 
mulation (Rutkowski et al, 2008), together with the fact that they 
were elicited using animals with a specific lesion in ER stress 
signaling, argue that they are likely to apply to ER stress in general. 

Two key questions immediately arise for further study. The 
first of these is whether all of the metabolic gene regulatory 
changes downstream of ER stress are subordinate to HNF4a 
and subsequent changes in the expression of just one or a 
small number of rapidly suppressed master regulators such as 



assessed by chromatin-IR Regions are given relative to the transcriptional 
start site, and correspond to regions identified by ChlP-seq analysis 
(Schmidt et al., 2010). n = 3-4 animals per group. Typical recovery of 
genomic material in samples containing HNF4a antibody was in the 
range of 0.1-1 percent of total input. *p < 0.05 by f-test. 
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FIGURE 8 | Model for organization of lipid metabolic gene regulation 
during ER stress. Experimentally demonstrated (here or elsewhere) 
relationships are shown using solid lines, while as-yet unvalidated 
relationships suggested by this work are shown by dashed lines. 



Cebpa, Pgcla, etc. Testing this hypothesis requires systematic 
overexpression of different metabolic transcriptional regulators, 
to see which blunt lipid accumulation during ER stress and how 
downstream metabolic genes respond. It is possible that HNF4a 
regulates both the proximal transcriptional regulators (C/EBPct, 
PGC-la, etc.) and also genes encoding the downstream metabolic 
enzymes (CPT1A, ACOX1, etc.), which would provide a feedfor- 
ward mechanism for suppression of these pathways. The second 
question is how ER stress influences HNF4a activity. Absent 
changes in HNF4a expression level or localization (Figure 7), the 
most likely scenarios are either covalent modification of HNF4a 
itself or modification, change in expression, or sequestration of a 
binding partner of HNF4a. 

It is possible that the preponderance of bZIP factors produced 
during ER stress (ATF6ct, ATF6f5, ATF4, XBP1, and CHOP) could 
titrate a coregulator away from HNF4ot. Persistent ER stress expe- 
rienced in Af/6a~/~ animals causes prolonged UPR activation 
and expression of the other non-ATF6a UPR-regulated bZIPs 
(Rutkowski et al., 2008), which might therefore cause continued 
cofactor sequestration even as wild-type animals recover from 
stress. A sequestration model predicts that the strength and per- 
sistence of the stress will be a key factor in altering the activity of 
non-UPR transcription factors like HNF4a. Accordingly, chronic 
stresses would be expected to elicit different patterns of metabolic 
gene regulation than acute stresses, and strong stresses different 
patterns from milder ones. In support of this idea, exposure of 
zebrafish larvae to chronic but mild TM results in efficient induc- 
tion of steatosis in the larval livers, and knockdown of ATF6a 
ameliorated this lipid accumulation. In contrast, a more acute but 
stronger ER stress — perhaps most akin to the exposures used in 
this work — led to less efficient steatosis that was exacerbated by 
ATF6a knockdown (Cinaroglu et al., 20 1 1 ) . The effects of stresses 
of different strengths and persistence have not yet been tested 
in a mammalian system, and will be important in validating or 
refuting the idea that pathophysiological conditions like obesity 



are in effect states of chronic ER stress. We also note that this 
work warrants an exploration of the roles of NR2F1 and ELK4 
in hepatic gene expression. In fact, with all genes broken down 
into 81 possible expression profiles, we predict that other testable 
hidden regulatory nodes will emerge, and that searching these 
groups for conserved sequences will reveal nodes beyond the rela- 
tively small number that are linked to transcription factors in the 
JASPAR CORE database. 

Finally, the scope of this work was limited to exploring the 
gene regulatory network linking acute ER stress to metabolism. 
Our gene expression results suggest that genetic suppression of 
VLDL production and fatty acid oxidation likely contribute to 
the steatotic phenotype while lipogenesis does not contribute, and 
might even be inhibited as a feedback mechanism to offset lipid 
accumulation. However, it remains to be determined whether 
these changes in mRNA expression are reflected in protein lev- 
els and consequent biochemical activities of the various pathways, 
and the extent to which they are also seen during physiological ER 
stresses such as obesity. 

In summary, we have provided proof-of-principle that a 
bottom-up approach sheds light on the organization of metabolic 
gene regulation during ER stress and makes testable predictions 
about this organization. As a consequence, we have identified a 
novel regulatory node in the process. Beyond revealing a likely 
role for HNF4a in this regulation, it provides a resource for 
regulators of other coordinated gene expression groups to be 
discovered. 
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